Dynamics of threads and polymers in turbulence: power-law distributions and 

synchronization 
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We study the behavior of threads and polymers in a turbulent flow. These objects have finite 
spatial extension, so the flow along them differs slightly. The corresponding drag forces produce a 
finite average stretching and the thread is stretched most of the time. Nevertheless, the probability 
of shrinking fluctuations is significant and is known to decay only as a power-law. We show that 
the exponent of the power law is a universal number independent of the statistics of the flow. For 
polymers the coil-stretch transition exists: the flow must have a sufficiently large Lyapunov exponent 
to overcome the elastic resistance and stretch the polymer from the coiled state it takes otherwise. 
The probability of shrinking from the stretched state above the transition again obeys a power law 
but with a non-universal exponent. We show that well above the transition the exponent becomes 
universal and derive the corresponding expression. Furthermore, we demonstrate synchronization: 
the end-to-end distances of threads or polymers above the transition are synchronized by the flow 
and become identical. Thus, the transition from Newtonian to non-Newtonian behavior in dilute 
polymer solutions can be seen as an ordering transition. 

PACS numbers: 47.27.Qb, 05.40.a 



I. INTRODUCTION 

Lagrangian chaos of the motion of fluid particles may 
lead to a non-trivial behavior of small objects dragged by 
the flow, because the exponential separation of nearby 
fluid particles, that holds in chaos, stretches different 
parts of small objects apart. While this is of no conse- 
quence if the immersed object is a solid, for elastic objects 
the stretching is important. Here we address the behav- 
ior of threads and of polymers immersed in a turbulent 
flow. Threads are defined as objects that resist the in- 
crease of their end-to-end distance only when the latter is 
close to the maximal length of the thread. For polymers, 
in contrast, the resistance exists at all distances, but it 
is only linear in a wide range of scales much smaller than 
the maximal length. The distribution of the end-to-end 
distances of these objects obeys a power law [lJ-Ql , while 
the exponent of the power law is non-universal Here 
we derive a new region of universality for the exponent. 
We also show that the flow synchronizes different parti- 
cles within spatial domains with extensions comparable 
to the Kolmogorov length, effectively making the vectors 
of the end-to-end distances to become equal. 

The main application of our results is to turbulence, 
where the positive Lyapunov exponent A is estimated as 
a characteristic value of the modulus of the velocity gra- 
dient Q and Lagrangian chaos holds. The details of the 
statistics of the turbulent velocity field are unknown, so 
one could expect that not much can be said on the behav- 
ior of the immersed threads and polymers. Nevertheless, 
we show that both a power-law behavior and synchro- 
nization hold. This is thanks to the universal nature of 
the derivation, which makes only quite general assump- 
tions on the statistics of the flow, that are plausible for 



turbulence. Our results are accessible to experiment and 
may be tested. 

Consider, placing a thread, with maximal length much 
smaller than the Kolmogorov length 77, in a turbulent 
fluid. Typically, the thread is coiled initially and the 
flow will straighten it by exponentially separating the 
thread's ends with time. The stretching will be arrested 
at about the maximal length. After that, the thread 
will randomly change its orientation in the flow with its 
end-to-end distance fluctuating near the maximal length. 
From the modern understanding of the behavior of poly- 
mers in the flow described below, it may be expected that 
the fluctuations are strong and that the frequency, with 
which small extensions occur in the flow, decays only as 
a power-law in the extension. Our results show that the 
exponent of the corresponding power-law in the distri- 
bution of the end-to-end distances is a universal num- 
ber independent of the statistics of the velocity (cf. @ 
and the discussion in the Conclusion). Thus we provide 
an analytic result for any turbulent (independent of the 
Reynolds number) or chaotic flow, that obeys the proper 
conditions of decay of correlations. 

A similar situation holds for polymers. For polymers 
the coil-stretch transition exists: when turbulence is not 
sufficiently vigorous, the elastic forces overcome poly- 
mer's stretching by turbulence and the polymer spends 
most of the time in the un-stretched, coil-like state. 
However when the strength of turbulent fluctuations in- 
creases, the stretching by the flow overcomes the elas- 
tic resistance and the polymer is stretched [Iol - [T2| . The 
quantitative criterium for the transition in the real tur- 
bulent flow was obtained in [l|. It was shown, with 
no modeling of turbulence, that the transition occurs at 
At = 1, where r is the polymer relaxation time, see also 
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Q. At At < 1, the polymer is coiled, while at Ar > 1 
it is stretched. A similar criterion was obtained indepen- 
dently for the white noise (Kraichnan) model of veloc- 
ity in 2]. It was also shown that the probability den- 
sity function of the polymer size obeys the power-law, 
both for real turbulence [l| and for the Kraichnan model 
For turbulence, the exponent of the power-law 
is non-universal, however near the coil-stretch transition 
it takes a universal form The numerical derivation 
of the power-law for turbulence above the transition was 
made in Q- However, the exponent could not be resolved 
numerically. 

Here we show for real turbulence, that the exponent 
of the power-law admits another universal limit. Namely 
we derive the universal expression for the exponent well 
above the transition at At ^ 1, cf. Q and the Comment 
below. Our prediction is accessible to the numerical test, 
cf. 0. 

An important property arises, when one passes from 
the consideration of a single particle to many particles. 
We show that the Lagrangian chaos orders the particles 
such that their end-to-end distances become identical. 
This gives a very important insight into the behavior of 
dilute solutions of such particles, allowing a better un- 
derstanding of the particles' back reaction on the flow. 
In particular, it gives some microscopic insight into the 
possibility to describe the polymer degrees of freedom 
above the coil-stretch transition by a single vector field, 
that was derived by macroscopic means in [l3l Il4| . The 
present work indicates that the respective vector, which 
gives the macroscopic equations of the dilute polymer so- 
lution a form similar to magnetohydrodynamics, is sim- 
ply the end-to-end vector of a single polymer, weighted by 
the number of polymers in the unit volume and the elas- 
tic properties of the polymer molecule. Furthermore, the 
fact that the effective form of the particles in the flow is 
rod-like, where "rods" of different particles are aligned, is 
important for deriving the correct criterium when the so- 
lution can be considered dilute. It was shown in that 
for aligned rods the condition of diluteness is much milder 
than for randomly oriented rods (the effective range of 
influence of one rod with another is ellipsoidal for the 
former and spherical for the latter). Here, we postpone 
the thorough discussion of implications of our results for 
the hydrodynamics of dilute polymer solutions to future 
work. We only notice that the implications are robust. 

These observations seem particularly important for 
turbulent dilute polymer solutions. More than 60 years 
ago, it was found that even an addition of a minute 
amount of polymers to a fluid can significantly change its 
properties such as turbulent drag ( "drag reduction" ) . 
Despite high practical importance of this phenomenon 
and much theoretical effort, the understanding of turbu- 
lence in dilute polymer solutions is still rather limited. 
It was clarified, though, that this is a threshold phe- 
nomenon: turbulence must be sufficiently strong to over- 
come the elasticity of individual polymers and stretch 
them. As the Reynolds number and A grow, for A larger 



than the coil-stretch transition value 1/t, the fluid be- 
comes non-Newtonian. What our result implies then is 
that this transition has similarity to a magnetic phase 
transition. The configurations of the polymers are uncor- 
rclated below the transition, while all polymers "point" 
in the same direction above the transition. Moreover, 
even below the transition there arise isolated ordered 
domains, within which the polymers are stretched and 
aligned. The volume fraction of domains, where poly- 
mers are stretched above a certain size, obeys a power 
law in that size. As At increases, these domains become 
more and more frequent in space, starting to occupy most 
of the volume as At crosses unity. 

II. SMALL THREADS IN CHAOTIC AND 
TURBULENT FLOWS 

We first consider the dynamics of a thread carried by 
a turbulent flow. We assume that the thread is much 
smaller than the Kolmogorov length rj, so the difference 
of velocities at the thread's ends can be approximated 
by the first term in the Taylor expansion. We study the 
dynamics of the end-to-end distance R, assuming that, 
as long as the thread is not stretched to about its full 
length Rmax, the ends are simply dragged by the flow 
with the local fluid velocity u. If x is the coordinate of 
one of the ends, then 

R = u(x + R, t) - u(x, t) « aR, R<Rm ax , (1) 

where <7y(i) is the matrix of velocity derivatives djUi 
taken in the frame of the fluid particle x(t). Without the 
constraint R -C Rmax, the above equation is the equation 
governing the separation of two Lagrangian trajectories 
(trajectories of fluid particles) in the flow. Hence R evolv- 
ing according to the above equation will become of the 
order of R m ax at some t. We assume that the thread's 
resistance to stretching can be described by a radial force 
and hence the equation on R can be written as 

R = u(x + R,t)-u(x,t)-VU(R), (2) 

with some "potential" U(R). We assume that VU(R) is 
appreciable only for R ~ Rmax-, where it grows indefi- 
nitely, reflecting the resistance of the thread to further 
stretching. For example, as a model for threads, one may 
consider the FENE model of polymers (l7j . 

eR 

VU F ene(R) = z T5y5 to, (3) 

1 — (H/Hmax) 

where the spring constant e is assumed to be small, such 
that the thread resistance to stretching is insignificant at 
R -C Rmax- Thus, a polymer with very large relaxation 
time (and, hence, well above the coil-stretch transition, 
see below) would be considered as a thread. 

Here we wish to describe the region R <C R ma x , where 
the statistics becomes largely independent of the details 
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of U(R). Decomposing R = R m ax exp[p]n, where h is a 
unit vector describing the thread's orientation, one has 

p = hah - f'(p), f'(p) = R^ ax e- p U' [R m ax exp[p]\ , 
h = ah — h(hah), (4) 

where we introduced the function f'(p) for further con- 
venience. For the FENE model we have 

Ifene{p) 



1 — exp[2p] ' 
!fene{p) = -| In (exp[-2p] - 1) . (5) 

It is essential that the dynamics of h decouples from 
the one of p, while p is driven by the "external", 
p— independent noise nan. Furthermore, the dynamics 
of h is independent of U(R), and is the same even with- 
out U(R), 



ah' — h (h ah), p' — h 'ah' . 



(0) 



This equation describes the evolution of the distance be- 
tween two infinitesimally-close trajectories of the fluid 
particles. The Lyapunov exponent A is the asymptotic 
logarithmic growth rate of the modulus of that distance, 



A 



lim - / h'ah'dt' 



(7) 



where the existence of the limit follows from the law of 
large numbers, provided a has a finite correlation time 
t c (which is the case for turbulence H). 

We conclude from the above that hah in Eq. (|4]) has a 
positive mean A. Then, at p <C — 1, where f'(p) is negli- 
gible, the dynamics is a random motion with a mean drift 
due to A. One has p — A + £(t), where (£) = 0. We des- 
ignate averages by angular brackets. From this equation 
one expects an exponential behavior of the steady-state 
probability density function (PDF) P ss of p, like the den- 
sity of a gas under gravity. We first study the problem 
within the so-called Kraichnan model, where the velocity 
field u is modeled as a Gaussian random field with zero 
mean and zero correlation time. Most of the results for 
the Kraichnan model can be inferred from , however 
there is one important new observation that we make 
here. The demands of statistical isotropy and homogene- 
ity in space imply that a is a stationary Gaussian matrix 
process with zero mean and pair-correlation function Q 

(a a p(t)a^ s (0))=D8(t) [{d+l)8 al 5p S -8 aS 8^-8 a p8^ s ] . 

Here, d is the space dimension. The above form respects 
the incompressibility condition and can be shown to be 
unique under the above assumptions on the statistics of 
u. For this model, hah is statistically equivalent |l8[ to 
the sum of A = d{d — 1)DJ2 and a white noise £, so the 
effective dynamics is (cf. [3|) 

P = \ + Z-f'(p), (Si(tW)) = (d-l)D6(t-t'). (8) 



Here, D is a diffusion coefficient. The associated Fokker- 
Planck equation gives in the steady state 

2 [(A - f'(p)) Pj - (d - 1)DP' S ' S = 0. (9) 

The steady-state solution must have a constant flux C, 

2[\-f'(p)]P ss -(d-l)DP' ss =C. (10) 

At p — > — oo, one can neglect the f'(p) term, and the 
solution becomes equal to a sum of C/A and an expo- 
nentially decaying solution. The condition for P ss to be 
normalizable gives C = 0, which leads to 



Pss(p) 



N 



1 ( 2Xp-2f{p) \ 
= — cxp I — — — I 



N 



exp 



(d-l)D J' 
2\p> - 2f(p>) 



(d-l)D 



dp', 



(11) 



where the properties of f(p) are such that the normaliza- 
tion integral iV is determined by R ~ R ma x or \p\ ~ 1- 
At p <C —1 one may neglect f(p), and the solution is 
exponential. At R <C R m ax, one finds 



Pss(R) 



NRr, 



R 



R n 



1 i 2A 
1+ (d-l)D 



R Rn 



It is tempting to consider X/D as a parameter in the 
above equation, hoping that this will allow to describe a 
wider class of physical situations Q. Nevertheless, it will 
be clear below that only the actual value of 2X/[(d— 1)D 
equal to d brings the physically meaningful result. After 
this substitution, the PDF for the end-to-end distance R 
becomes 



Pss(R) 



1 



R 



d~i 



R Rn 



A^ Rmax \ Rmc 

In particular, for the FENE model one finds 



.(12) 



P. 



FENE (p) = lc Xp (pd-& 



x (1 - exp[2p]) 
For P SS (R) we find 

1 



N 

de/[2\] 



(13) 



P. 



FENE 



(R) 



1 - 



R 



NR. 

2 



1? 



Rr< 



max \ Rmax 
-i de/[2\] 



d—l — pde/ A 



(14) 



We observe in Eq. (fT2"j) that the PDF is a positive power 
law (the model assumes incompressibiliy and is mean- 
ingful only for d > 2). Thus, in the limit of zero corre- 
lation time, the PDF of the end-to-end distance decays 
to small R rather slowly. Furthermore, the PDF pT|) 
depends only weakly on the details of the non-elasticity 
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expressed by f(p) and of the strength D of the velocity 
gradient fluctuations via a multiplicative constant. For 
larger D, the flow stretches the thread stronger, but the 
fluctuations of hah to negative values are also stronger. 

We test the above predictions numerically for two and 
three dimensions, using the following discretization (cf. 
[20l| ^ . We employ the random renewal model, where 
<?ap(t) is a piecewise continuous process, a(t) = a k at 
kAt < t < (k + l)Ai, and where a k are random, in- 
dependent Gaussian matrices with zero mean and pair 
correlation function 



Time fraction at a given extension vs. theory 



6 kl D 
At 



[(d + l)S a7 S/3s — S a 55p 7 — 5 a p5 7 s] 



The above process reduces to Eq. (Q} in the limit eo = 
DAt — > 0. For the three-dimensional case, the Gaussian 
matrix with the above correlation function and zero mean 
can be generated using the ansatz 



>afl 



1 

At 



VZe aM V k + V3f k f k - 



k\2 



V3 



where e Q /3 7 is the antisymmetric symbol, and V k and 
f k are independent Gaussian random vectors with zero 
mean and pair correlations (V k Vg) = 8 kl 8 a p and 

(/a/a) — S kl 5 a /3, respectively. In the two-dimensional 
case, the ansatz is 



'a/3 



1 

A* 



2e afS3 V k 



V2f k f$ 



k\21 



V2 



where V k is a set of independent random variables with 
zero mean and pair correlation (V k V l ) — 5 kl , and f k are 
independent two-dimensional Gaussian random vectors 
with zero mean and pair correlation (f k f l p) = 8 kl 8 a p. 
The above formulas make it clear that in the Kraichnan 
model the antisymmetric part of er, which describes the 
vorticity, and the symmetric part, which describes the 
strain, are independent. We also note that in two dimen- 
sions the evolution matrix exp[Ata k ] could be written 
explicitly, by rewriting the ansatz above with the help of 
Pauli matrices <7j, 



Ata k = ^hhai 



/I) *s/2. 



After V k (V k ) and f k are generated, the evolution of the 
end-to-end vector R k — R k h k of the thread follows from 



xk+l 



N k+1 = exp[Ata k ]h 
R k+1 = mm[R k N k+1 ,R Tnax ], 



jV fc + 1 



(15) 



where R k N k+1 = \exp[Ata k ]R k \. The last equation of 
([TT5j) describes the simplest model of resistance to stretch- 
ing: the thread is stretched as if its ends are two fluid 
particles, unless the distance between the ends reaches 
Rmaxi a t which value the distance growth is arrested. 
This corresponds to a potential U (R) , which rises sharply 
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FIG. 1: Steady-state distribution for the end-to-end distance 
R of a thread immersed in a turbulent fluid in three dimen- 
sions. R on the abscissa is given in multiples of the bin width, 
1/1200. The distribution is not normalized. The smooth line 
is a parabolic fit, as predicted by Eq. (|12[) for d — 3, to the 
experimental points. 



near R m ax- Hence, one may expect that the power law 
will extend up to values close to R ma x- 

The numerical result of 10 5 iterations of the above dy- 
namics in three dimensions (where, without loss of gen- 
erality, we use reduced units, for which R m ax = 1 and 
eo = DAt = 0.02) is shown in the histogram of Fig- 
ure [TJ There, the interval < R < 1 is split into 1200 
equally-spaced bins labeled by i, and the number of oc- 
currences of R k in each bin is counted. This sum is pro- 
portional to the fraction of time that R spends in the 
respective bin, for which the theory, Eq. (Ti~2"l) . predicts 
a parabolic dependence on R, for R sufficiently smaller 



than R r , 



1. Removing 50 b ins near R — R n 



1, the resulting histogram shown in Figure Q] fits the 
parabola Ci 2 ,i G {1, 2, . . . , 1150} for C = 1.2 x 10~ 4 
very well (the scattered points deviate from the fit by 
less than 10 %). The power law even extends to values 
of R very close to R m ax in accord with the expectation 
mentioned above. Analogous results for simulations in 
two dimensions are shown in the histogram of Figure [21 
where 2 x 10 5 iterations were made. A very satisfactory 
fit to a linear profile, as predicted by Eq. (TT2"j) for d = 2, 
is shown by the straight line, where the scatter of the 
experimental points does not exceed 10 %. The physical 
significance of the power law is illustrated in Figure [3] for 
the three-dimensional case. It shows a realization of R k 
in the steady state over 1000 steps of its evolution. It 
is clear that the excursions to small R occur quite often, 
as the slow power-law decay of P ss (R) to small R would 
make one to expect. 

It is possible to generalize the above results to any 
finite-correlated er(t), in particular, to the one of turbu- 
lence. Such a generalization is very important, as it gives 
an experimentally testable prediction for the behavior of 
threads in a chaotic flow and, in particular, in turbu- 
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Time fraction at a given extension vs. theory 
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FIG. 2: Steady-state distribution for the end-to-end distance 
R of a thread immersed in a turbulent fluid in two dimen- 
sions. R on the abscissa is given in multiples of the bin width, 
1/1200. The distribution is not normalized. The smooth line 
is a linear fit, as predicted by Eq. (|12[1 for d = 2, to the 
experimental points. 



Realization of the end-to-end distance in the steady state 
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FIG. 3: Shown is a time series R for 1000 time steps of a 
thread in three dimensions. The excursions to R significantly 
smaller than R m ax = 1 are quite frequent. 



lence. It requires to take the Lyapunov instability of the 
flow into account, which involves the time evolution of 
infinitesimally-close trajectories. This is considered in 
the following section. 



III. SEPARATION VECTOR BETWEEN TWO 
INFINITESIMALLY-CLOSE TRAJECTORIES IN 
A CHAOTIC FLOW: A NEW RESULT 

Below we derive two results on the behavior of the sep- 
aration vector between two infinitesimally-close trajecto- 
ries in a chaotic flow. The first result is new, while the 
second one is less known. These results are at the basis 
of several physical predictions obtained in the following 



sections. Where there is no ambiguity, we use the same 
notation as used in other sections for other quantities. 



A. Relaxation of the orientation of different 
separation vectors to the same vector 

It follows from the previous sections that the distance 
R between two infinitesimally close trajectories obeys 



dR 



R(t) = W(t)R(0), 



(16) 



where we introduced the Jacobi matrix W(t) of deriva- 
tives of fluid particles' coordinates with respect to their 
initial position. If q(t, x) is a Lagrangian trajectory, 



dq{t,x) 
dt 



u[t,q(t,x)], q(0,x) = x, (17) 



then the Jacobi matrix describes the distance R(t) be- 
tween two infinitesimally close trajectories starting at x 
and x + R(Q) according to 



R(t) 

W l3 {t,x) 



q(t,x + R)- 
dqi(t,x) 



q(t,x) « W(t,x)R, 



(18) 



Below we suppress the spatial index in W(t, x), where it 
is clear from the context. The matrix W(t) obeys 



^- = ° W > W(0) = 1, 



(19) 



in accord with Eq. (TTBl . The variable p(t) and the unit 
vector n(t) defined by R(t) = R(0) exp[p]h(t) satisfy 



dn 
~dt ' 

p(t) 



an 



handt , 



h(t) 



lim — 

t— >oo t 



W(t)h(0) 
nahdt = X. 



(20) 



We introduce the matrix decomposition W(t) = RAN, 
where R and N are orthogonal matrices, while A is diagi- 
nal, Ajj = Sij exp[Xi(t)t]. According to the Oseledec the- 



orem [21|, lim^oo Xi(t) = A,, where Ai are the Lyapunov 
exponents, while N(t), diagonalizing W T W, saturates at 
large times at a constant matrix N°°. Here, the index 
T means transposition. We assume that the Lyapunov 
exponents are ordered, A^(t) > A i+1 (£), and that Ai is 
strictly greater than A2. We have 

[W{t)h{Q)]i = Rtj exp[X 3 it)t]N jk n k {0), (21) 

where the summation over repeated indices is presumed. 
At times, such that exp[Ai(i)i] ^> exp[A2 (£)£], the j = 1 
term dominates the sum over j. Introducing the unit 
vectors rhi = Rn and U = Nn, we find that 



W{t)n(0) « exp[Ai (t)t}m(t) I ■ n(0) 
W{t)h(0) 



h(t) 



\W(t)h(0)\ 



m{t)sign I ■ n(0) , (22) 
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where the sign factor corresponds to interchanging the 
labels of the trajectories and has no physical significance 
for the applications we consider in the following sections. 
The approximation above holds for all «(0), for which I ■ 
ri(0) is much larger than the small parameter exp([Ai(i) — 
\2(t)]t). The Oseledets theorem implies that I saturates 
at a constant vector Iq at large times [2l|. This can be 
verified directly by deriving from the equation for W the 
following equation for 1, see 



d l t J^ (R T crR) lj + (R T aR) 

rH 2.^ s i n V, \f(\, (i\ - Xj(t))] 



./I 



dt ^ sinh[£(Ai(t) 

3 — * 



Nn. (23) 



At times considered in Eq. (|2"2")l , the time derivative 
of I decays exponentially, which indicates that I is ap- 
proximately constant at those times. Thus the approx- 
imation (|2"2"|) holds for all initial conditions besides the 
exp[(A2 — Ai)i] vicinity of the plane Iq ■ n(0) = 0. 

Choosing with no loss the sign of n(0) so that l o -n(0) > 
0, we find from Eq. that 

n(t) m(t), cxp[(A 2 - Xi)t] « min[l, l ■ ra(0)].(24) 

It is possible to verify from the equation for R(t) that 
m(t) obeys the same equation as n(i), so the above ap- 
proximation is consistent [l4j |. Since the corrections to 
Eq. (pM)) are of order exp[< (Ai — A 2 )], we conclude that 
different n(t), obtained from different initial conditions, 
exponentially relax to the same unit vector to. 

The most remarkable feature of Eq. (|24)) is that n(t) 
lost all dependence on the initial condition n(0), provided 
?i(0) is not in the narrow vicinity of the plane lo-n(Q) = 0. 
Thus, any initial condition outside the latter plane will 
relax to rh(t) at sufficiently large times. 

The observation of this property seems new, and we 
consider it now in more detail. We introduce a unit vector 
q(t) by q(t) = N(t)h(0). The Oseledets theorem ensures 
that q(t) saturates at a constant vector q°° = 7V°°fi(0) 
in the limit t — > oo. The following exact relation follows 
from Eq. ([20]), 



R T (t)n(t) = 



mm 

\Aq(t)\ 



exp[t(A 2 (i)-Ai (*))], 



■exp^A^-A^i))] , 



(25) 



" ' " ' Q 

where we defined 

Q 2 = q\ + q\ cxp[2t (\ 2 (t) - Ai(t))] + . . . 
+ exp[2t(X d (t)-X 1 (t))]q 2 d . 

We find that for such n(0), for which q^° ^ 0, we have 

lim R T (t)h(t) = (1,0, ...,0), (26) 

t— »oo 

where it is assumed that the initial choice of the labels 
of the trajectories corresponds to the positive sign of Iq ■ 



n(0) (otherwise there appears —1 in the vector on the 
RHS above instead of 1). Since m(t) = R(l, 0, . . . , 0), 
the above equation gives the more precise formulation of 
Eq. (pM)) . The condition qf ^ means that the limit 
holds for all h(t) such that n(0) obeys Iq ■ n(0) ^ 0. This 
limit holds for almost every trajectory, with the possible 
exception of a set of trajectories with a vanishing volume 
of their initial conditions. 



B. A relation for the exponential growth of the 
distance between two innnitesimally-close 
trajectories 

It follows from the Oseledets theorem that the sepa- 
ration R(t) between two infinitesimally close trajectories 
behaves exponentially, 



lim I In (44) = A. 



(27) 



For the case of incompressible flow considered here, the 
above limit holds for any initial position of trajectories, 
possibly except for a set of positions with zero spatial 
volume. Thus, one could expect R a (t) to behave as 
exp[aAi]. This expectation fails due to the very inter- 
mittent nature of the growth of the moments. There 
are small regions in space, for which R(t) grows with a 
time-dependent (local) exponent larger than A, and other 
regions, where R(t) grows with a time-dependent local 
exponent smaller than A, or even decays exponentially. 
The volume of these regions exponentially decays to zero 
as t — > oo, and it can be described as follows. One in- 
troduces x(t) = In ("j^gy) A an d describes it by the PDF 

P(x,t), where the probability is defined by the fraction 
of the volume for which x(t) = x. Then the exponential 
decay of the volumes, for which x(t) differs from A, is 
described by 



P(x,t) ~ exp[-tS(x - A)] , 



(28) 



where S (x) is the rate function, which is convex and pos- 
itive everywhere except at zero, see Ref. Q and the ref- 
erences therein. Although for large times the volume of 
the regions with x ^ X is exponentially small, these re- 
gions still contribute significantly to the spatial averages 
of R a (t). This is, because for these regions R a (t) can 
be exponentially larger than exp[aAi] (or smaller as may 
be needed for a < 0). The resulting non-trivial statis- 
tics is described by the exponential growth function 7(a) 
defined by @,[22[, 

7(a) = lim - \n(R a (t)) = lim - ln(\W (t)h\ a ) , (29) 

t— >oo t t— >oo t 

where n — n(0). Using the cumulant expansion, one can 
show that 7(a) is well-defined, and using P(x,t), one 
can see that 7(a) is the Legendre transform of the rate 
function S(x). The Lyapunov exponent enters 7(a) via 
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the relation 7'(0) = A, as can be easily inferred from 
the definition. Furthermore, it follows from the Holder 
inequality that 7(a) is convex (which is equivalent to the 
convexity of S(x)). 

Restoring the spatial argument of W{t), see Eq. (fl"8|). 
one can write down the average explicitly, 



(\W(t)nn 



^\{h-V m )q(t,x)\ a , (30) 



where O is the system volume. For isotropic statistics the 
average above is independent of h. Since the statistics of 
a is determined by the small-scale turbulence, where the 
statistics gets isotropic, we assume below that isotropy 
holds. 

We consider the integral of |Wn| over the direc- 
tions h and make the following change of variables in 
the expression for the average, n' = Wh/\Wn\. This 
transformation was introduced in Ref. (22j , where it was 
shown that the Jacobian of this transformation of the 
unit sphere is |tUn| _d , i. e. dn' — \Wn\~ d dn. We find 



dh 



\W(t)h 



= / \W-\t)h'\»dh', (31) 



where we note that IIU" 1 ?^'! = l/|IUri|. Averaging 
that relation, which holds at every point [remember that 
Wij(t,x) = djqi(t,x)], and using the fact that the aver- 
age is independent of the direction, we obtain 



(\W(t)n 



(\W-\t)n 



(32) 



zeros are given by 7'(0) = Ai and j'(—d) — Ad- These 
are positive and negative, respectively, provided the sys- 
tem is non-degenerate and Ai > (then incomprcssibility 
implies that A^ < due to J^A^ = 0). It follows then 
from the convexity that 7(a) < for — d < a < 0, and 
7(a) > for a > and a < —d. 



C. The exponential growth rate function for the 
Kraichnan model 

The expression for 7(a) within the Kraichnan model is 
well-known, see [8j and references therein. Here, to keep 
the exposition self-contained, we rederive the result. Us- 
ing the representation R(t) — R(0) exp[p(t)]n(i) , where 
\h\ = 1, we find that the definition (|29|) of 7(a) admits 
the following form: 



7(a) = lim — In ( exp 

t->-oo t 



a 



h(t')a(t')h(t')dt' 



(36) 



For the Kraichnan model n(t')a(t')n(t') is Gaussian, and 
it follows that 7(a) is parabolic. From j(—d) = 7(0) = 
and 7'(0) = Ai = A one infers that 



7(a) 



\a(a + d) 



(37) 



Note that for the Kraichnan model A^ = j'(—d) = —X\. 



The above relation allows to draw an important con- 
clusion for 7(a). First, putting p = 0, we recover the 
relation (\W(t)h\~ d ) = 1, that was established in [22j]. 
Second, differentiating the above relation with respect to 
p and substituting p — in the resulting relation, we 
obtain 

(\W(t)n\- d \n\W(t)n\) = -(In |W _1 (*)^D- (33) 

Using the definition (|2l?]) of 7(a) to write down 7 / (— d), 
and then using the relations above, we obtain 



l'(-d) 



lim 

t— >QC t 



1 (\W(t)n\- d ln\W(t)h\) 



(\W(t)r 



lim -(lnlFy-^nl). 

t— >oo t 



(34) 



However, the last term is nothing but A^, as the aver- 
age of the logarithm is determined by \d(t) due to the 
appearance of W~ l in the average. We find 



j'(-d) = A d 



(35) 



This relation can be found in the Appendix of [14]. It is 
less known but important for our analysis below. 

We now summarize the properties of 7(a), that will 
be important in the following. This function has two ze- 
ros, 7(0) = 7(— d) = 0. The derivatives of 7(a) at these 



IV. DERIVATION OF P SB {R) oc i? d 
TURBULENCE 



FOR REAL 



We now generalize the result on the power-law distri- 
bution of R at R <C Rmax to any finite-correlated cr(t), 
in particular, to the one of turbulence. This gives an 
experimentally testable prediction on the distribution of 
the end-to-end distance of a small thread in turbulence. 
We then use the technique introduced in [l| for a study 
of the distribution of polymers in turbulence. 

At p <^ —1, over a not too large time interval t such 
that p(r) <C — 1 for < r < t, we have 



p(t) = p(0) + z(t), z(t) 



nafidt' . 



(38) 



Generally speaking, z(t) grows with time, and at a suffi- 
ciently large t one can neglect the contribution to z(t) of 
a r c neighborhood of t = 0. Then z(t) is approximately 
independent of p(0) that, by causality, is determined only 
by a(t) with t < 0. It follows that the PDF of p(0) + z(t) 
is given by the convolution of PDFs of p(0) and z(t). In 
the steady state both the distribution of p{t) and of p(0) 
are given by P ss (p). It follows that P ss (p) must obey the 
following stationarity condition, 



Pss(p) = (Pss [P-Z{t)\). 



(39) 
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Using Laplace transforms, one can show that the above 
equation has exponential solutions oc exp[/ip], where p 
obeys (cf. fl) 



exp 



-fi I nandt' 
'o 



= 1. 



(40) 



It follows from the analysis in the previous section that p 
obeys 7(— p) = 0. There are two solutions, which are the 
same as in the Kraichnan model: /i = corresponding 
to P ss (p) — const and p — d corresponding to P ss oc 
exp[/5<i]. Only the latter solution is normalizable at p — > 
— oo. Passing to the variable R, we obtain 



Pss(R) 



1 



NR n 



R 



Rr, 



d-l 



R <c R n 



(41) 



where N is a constant determined by matching the so- 
lution above to its form at R ~ R ma x- The latter can- 
not be obtained for a general a. We stress that although 
Eq. (|41l) gives the same result as in the Kraichnan model, 
this does not signify that an effective white noise descrip- 
tion is possible in the region R <C R ma x, as it will become 
more clear in the next sections. Furthermore, the Kraich- 
nan model gives the correct result only for the ratio X/D 
that holds for the full model. This is most essential for 
the understanding of the applicability of the model, cf. 

Thus, the dynamics of a small thread thrown into a 
turbulent flow consists of stretching to its full length 
and rotation as a rigid rod, which is intermittently inter- 
rupted by the shrinking and relaxation of the thread. The 
shrinking intervals - excursions to the domain R < R ma x 
- occur rather often in accord with the power-law distri- 
bution proportional to R d ~ 1 . To check the theoretical 
prediction, one needs to measure the PDF of a small 
thread in a three-dimensional turbulent (or chaotic) flow 
and compare it with an .Redistribution, or make the mea- 
surement in two-dimensions, and compare the result with 
R. 



V. SYNCHRONIZATION OF THREADS BY 
THE FLOW 

Having described the distribution of the end-to-end 
distance of a single thread, let us compare the behavior 
of different threads in the same velocity gradient a. For 
turbulence this means that the threads are separated by 
a distance much smaller than the Kolmogorov length m 
but not too small, so that the hydrodynamic interactions 
between the threads are still negligible. We consider the 
dynamics of the orientation vector n. As the dynamics 
of n is the same as for two fluid particles, see Eq. fp~6|) . 
we find immediately from Eq. ([24)) that orientations of 
different threads at large times are equal to the same 
vector, 



n(t) ~ m(t). 



(42) 



Here, we omit the sign factor [see Eq. (|24"|) ]. as a change 
in the sign of h does not change the configuration of the 
physical thread in space. Thus, threads with initially 
different orientations relax to the same vector exponen- 
tially. As it is clear from our considerations in Section 
IIII Al the exponent of the relaxation can be estimated as 
Ai — A2. 

One can expect that not only the orientations, but also 
the extensions of the threads relax to the same value, such 
that the flow synchronizes the threads completely, and 
the end-to-end vectors become equal for different threads 
in the same flow. Consider the radial variables p\ and 
P2, each of which obeys the first of Eqs. ((4]). Due to the 
equalization of orientations, at sufficiently large times the 
dynamics reduces to 



= mam - f'(pi), 



(43) 



so that - if we introduce bp 
of two equations, 



: P2 — Pi - we obtain a system 



dSp ,,, , 
-& = f (Pl) 



f'(px+Sp), 



Pi 



f'(pi) 



This system may be considered as a Langevin-like sys- 
tem of equations, with an obvious-steady state solution 
P ss (Sp, pi) — 5 (Sp) P ss (pi), where the steady-state solu- 
tion for a single radial variable, P S8 (p), was discussed be- 
fore. One expects that this solution provides the unique 
steady-state solution for the system, signifying complete 
synchronization between different threads. For example, 
for the FENE model one has 



dSp 



1 



1 



1 — exp[2pi] 1 — exp [2pi + 26 p] 



The time derivative on the RHS always has the sign oppo- 
site to that of dp (remember that both p\ and pi — pi+Sp 
are negative) . The above equation describes monotonous 
relaxation of dp to 0, which demonstrates synchronization 
for the FENE model. We also verified the synchroniza- 
tion numerically for the three-dimensional renewal model 
of Sec. |TTJ We constructed the Jacobi matrix for n evo- 
lution steps, 



W n 



n° xp 

fc=l 



rk rk 



6<W(/ fc ) 2 



and verified that 



[W n WZ]ij RiiRji 



fiihj exp[/?i 



TliUj 



(44) 



relaxes to a matrix composed only of ones for different 
initial conditions of h. The relaxation of all elements 
to unity (up to order 10~ 4 ) was already observed after 
50 — 100 evolution steps, slightly depending on the ini- 
tial conditions. This verifies the alignment of different 
threads as described by Eq. (|4"2"j) . Then we verified that, 
with the same number of evolution steps, the values of 
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the radial component p relax to the same value for dif- 
ferent initial conditions. Since for the Kraichnan model 
for d = 3 one has A2 = and Ai — A2 = A = 3D, see 
Ref. [8]. which in our dimensionless units equals 3e with 
e = 0.02, the convergence occurs within the expected 
time-scale (Ai — A2) _1 . It seems highly unlikely that 
our conclusion on the relaxation of the radial coordinate 
is model-dependent, and we will assume that it holds 
for any physically meaningful f(p). Thus, our analysis 
brings the conclusion that different threads get synchro- 
nized by the flow. This result can be tested experimen- 
tally. 



VI. THE EXTENSION DISTRIBUTION FOR A 
SINGLE POLYMER ABOVE THE 
COIL-STRETCH TRANSITION 



As mentioned above, threads can be considered as the 
limit of a polymer with a very large relaxation time, 
which is the case of the FENE model for the thread. One 
may expect that the results are generalizable to polymers 
with an arbitrary relaxation time, which we pass to show. 
The main difference of polymers from threads is that the 
former exert a significant elastic resistance force, which 
opposes stretching already for R <C Rmax- In particular, 
this force may arrest the stretching by the flow, which 
corresponds to the regime below the coil-stretch transi- 
tion [l|, [nj-Qjf • The contents of this section and the next 
one are largely inferrable from the previous work and are 
introduced to keep the exposition self-contained. 

The results for threads can be carried over to polymers 
above the coil-stretch transition rather straightforwardly. 
The following effective equation holds for the end-to-end 
vector R of the polymer (see e. g. Ref. |2[), 

R = aR- — -VU(R) + £, (45) 

T 

where r is the polymer relaxation time, and V£7 (R) is 
negligible for a wide range of scales R co u <C R <C Rmax- 
The white noise £ produces equilibrium fluctuations, that 
govern the behavior of the polymer in the coiled state 
where R ~ R co u "C Rmax- This noise can be neglected 
for Ar > 1 above the coil-stretch transition (see Ref. 
[l|), where the polymer is stretched to R >• R co u and 
£ is small in comparison with the term R/t. The self- 
consistency of this omission is checked by dropping £ in 
the equation above and passing again to the same vari- 
ables p and h introduced before for threads (with no am- 
biguity we use here the same notation). One finds 

p = hah f'(p)> ^ = an ~ h[hah). (46) 

T 

Like for threads, the dynamics of the orientations is the 
same as for a vector separating two fluid particles. Thus, 
the only change in comparison with threads is that the 
Lyapunov exponent should be substituted by the effec- 
tive value A — r -1 . For Ar < 1, this change is essential, 



as the above equation - on average - describes a p that 
decreases for R ^> R co u, giving inconsistency if C is ne- 
glected. Above the coil-stretch transition, however, at 
Ar > 1, the neglect is justified. In the Kraichnan model 
0-111], the straightforward generalization of the analysis 
for threads yields for R <C R max 



Pss{R) 



1 



N'Rr, 







( J , a-d 









,(47) 



where N' is a constant, which should be determined by 
matching the above solution to the asymptotic region 
R ~ Rmax- The above result is the same as for the FENE 
model of the thread with e = 1/r. This is necessary, 
because the calculation for the threads did not use the 
smallness of e. Below the coil-stretch transition, a < 
and the above PDF is not normalizable at small R, which 
indicates the inconsistency of neglecting £. The power 
law below the transition was derived in Refs. [H, 0. In 
contrast, above the coil-stretch transition, a > and the 
distribution is normalizable at small R, indicating that 
the normalization is determined by R ~ Rmax, where the 
polymer spends most of the time. Since the derivation 
of the power law is based on considerations local in R, it 
is readily concluded that the power law holds at R co u <C 

R Rmax- 

A noticeable difference from the case of the thread is 
that for polymers the exponent of the power law is non- 
universal and depends on the parameter D of the Kraich- 
nan model. This is due to the presence of the dimension- 
less parameter Ar. It should also be noticed that even 
above the coil-stretch transition, P ss (R) can be a decreas- 
ing function of R, which has a normalizable singularity 
at R = 0. 



VII. THE GENERALIZATION TO THE CASE 
OF REAL TURBULENCE 

The generalization to the arbitrary statistics of a, in 
particular to the one of turbulence, is a straightforward 
repetition of the analysis made previously for the threads 
and it follows [l| , see also [23| . One has at p <c — 1 that 



p(t)=p(0) + z(t)--, z(t) 

T 



hahdt'. (48) 



Again, for sufficiently large t, one neglects the contri- 
bution to z(t) of a r c neighborhood of t = and finds 
the following condition for the steady-state distribution 

Pss(p)- 



Pss(p) = P, 



p 



z{t) + - 

T 



(49) 



Again, this equation has exponential solutions oc expl/tp], 
where p obeys (cf. [l[) 



exp 



— \i \ hahdt' 


\ = exp 


pt 


Jo 




T 



(50) 
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The exponent p, obeys in this case the condition 
7(-/x) = -h/t. 



(51) 



As before, the convexity of 7(a) (and thus of j(— p,) + 
p/r) implies that there are just two solutions to this con- 
dition. Also again, the solution p = should be dis- 
carded, as it is non-normalizable. The exponent p, how- 
ever, is not universal such as in the case of the threads, 
because the behavior of 7(a) in general is not universal, 
but depends on the details of the statistics. The only gen- 
eral implication possible is that p obeys p < d, which fol- 
lows from the fact that 7(a) is negative for — d < a < 0, 
and positive otherwise, cf. Ref. [22j ]. 

Thus, the steady state of the polymer molecule in the 
flow above the coil-stretch transition satisfies 



P SS (R) 



1 



N"R. n 



R 



Rn 



(52) 



where 6 = 1 — p/d > 0. For the Kraichnan model 
5 = 1/(At). The condition of self-consistency of the 
above equation, guaranteeing that it describes polymers 
above the coil-stretch transition, is <5 < 1 which, for the 
Kraichnan model, reproduces the universal criterium for 
the coil-stretch transition, Air = 1 (see Ref. [l[ and 
below) . 



VIII. FROM POLYMERS TO THREADS: THE 
CASE OF A VERY LOOSE POLYMER 

One case, for which the general analysis provides the 
exponent governing the steady state distribution P ss (p), 
is the thread limit At ^ 1. We have in this case that p 
is close to the solution for the thread 



p = d[l-S], 5<1. 



(53) 



We can find 5 to leading order in 1/Ar, by using the 
relation 7 '(— ci) = Xd- The expansion of 7 (— p) near p = d 
in the steady-state condition 7(— p) = —(a/t gives 

j(-d + dS) ps d\ d S --, S |A d |r « 1. 

T A d T 

Thus, well above the transition, the exponent p depends 
linearly on the inverse relaxation time (remember that 
X d < 0), 



p = d+^-, |A d r|>l. 



(54) 



Correspondingly, the steady state distribution of the 
polymer extension well above the transition satisfies 



Pss{R) 



1 



N"R n 



R 



Rn 



dll-IAdrr 1 ]-! 



(55) 



For the Kraichnan model |A<j| = A, and the result from 
above holds generally. This prediction gives the power- 
law distribution for polymers in a flow well above the coil- 
stretch transition. Both terms in the above expansion 



constitute a new result that can be tested experimentally, 
either by putting polymers with a large relaxation time 
in a given flow, or by considering a vigorous flow with 
high I Ad I for a polymer with a given r. 



IX. THE DISTRIBUTION IN THE VICINITY 
OF THE COIL-STRETCH TRANSITION 

Another limiting case, opposite to the situation well- 
above the transition just described in the previous sec- 
tion, concerns the vicinity of the coil-stretch transition, 
for which 1/r pa Ai. The corresponding result below the 
transition was found in [l|), and while the generaliza- 
tion to the case above the transition is straightforward, 
we bring it here for completeness. It is clear from the 
stationarity condition (1511) that p is close to zero in this 
case. Expanding j(—p,) near p, = to second order, we 
find 



-\ x p 



7"(0)M 2 



T 



Discarding the trivial solution /i = 0, one has 

2Ai 2 7 '(0) 



p 



1 



1 



7"(0) 7"(0) 



(56) 



(57) 



It is clear from 7 (— d) = that k, which is positive due to 
the convexity of 7(0;), is estimated as d. This agrees with 
the result of Ref. [l[ , where it was used slightly below the 
coil-stretch transition. Here we only show that the same 
result also applies slightly above this transition. We have 



fi = K-— , |l-Air|<l. 

AlT 



(58) 



This clearly shows that the coil-stretch transition occurs 
at Air = 1, where fi changes sign. The corresponding 
steady-state distribution of the polymer extension near 
the transition satisfies 



Pss(R) 



1 



N"R ri 



R 



Rn 



Kll-i^r)- 1 ]-! 



(59) 



For the Kraichnan model, we have 7"(0) = 2Ai/d, see 
Eq. (f3"Tl) . Thus, k — d and fj, = d[l — (At) -1 ], in agree- 
ment with the previous results. The numerical confirma- 
tion of the result can be found in [24l n25| . 



X. SYNCHRONIZATION OF POLYMERS IN 
THE FLOW 

The synchronization of polymers in the flow is con- 
sidered along the same lines as for threads. The mod- 
ification of the motion equation for the difference Sp of 
the radial variables in Sec. |V]only concerns the addition 
of the term —Sp/r, which makes the relaxation of the 
extensions to the same size even faster. 
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There is, however, one important new regime, that ex- 
ists for polymers, but not for threads: the regime below 
the coil-stretch transition, where At < 1. There, the vol- 
ume fraction of the regions, where polymers are stretched 
to a size not smaller than some R satisfying R 3> R co u , is 
much smaller than unity and behaves as a positive power 
law in R co u/R. Still, these regions are rather frequent 
in space, as their decay is only a power law. The theory 
constructed here then shows that the polymers are syn- 
chronized in these rare regions, where they are strongly 
stretched. It was shown in Ref. [l[ that the rare fluctua- 
tions of the velocity gradient a, that form R much larger 
than the most probable size R co a, persist during a long 
time-interval, where one can apply the dynamics 



dR 



R 



(60) 



This dynamics holds during a time-interval much larger 
than (Ai — A2) _1 . Thus, these rare fluctuations of the 
flow, which stretch the polymers, also synchronize these 
polymers by starting from initially uncorrelated coiled-up 
states. 



XI. SYNCHRONIZATION AND 
HYDRODYNAMIC EQUATIONS OF DILUTE 
POLYMER SOLUTIONS 

The main application of the phenomenon of synchro- 
nization is to the hydrodynamic equations of dilute poly- 
mer solutions. The synchronization implies that in the 
regions where the polymers are stretched, one can intro- 
duce a smooth macroscopic field R(x,t), providing the 
end-to-end distances of polymers near the point x. This 
field is defined by 



R(x,t) = ^ f R i (t)5[x-x i (t)], 



(61) 



where i is the index of the polymer, and Xi(t) is the 
position of the molecule's center of mass. Differentiation 
yields 

-£-+u- VR = R-Vu~ — , (62) 

Ot T 

where one writes the substantial derivative of R(x, t) in- 
stead of the ordinary time derivative appearing in the 
equation for J? of a single polymer. As it was observed 
in [lij , the vector field equation of the above form leads to 
the purely decaying equation (9t+u-V)V-i? = — V-R/t 
for V-JfJ. Thus, in the steady state one can always assume 
V • R — 0. This conclusion on the spatial distribution of 
polymers seems to be non-trivial. 

The above equations imply that, if the polymer density 
n obeys the equation (we use V • u = 0) 



On 
~dt 



u ■ Vn = KV 2 n, 



(63) 



where k describes the diffusion of polymers, then the field 
Iljj = nRiRj obeys 

^ _|_ u . vn = iivm + (VtiVn - — + KR l R j v 2 n. 

dt t 

In many situations one can neglect the last term. In 
particular, this is true when polymers are uniformly dis- 
tributed in space, which is a possible steady-state solu- 
tion. It should be mentioned that there are indications 
that in certain experimental situations the inhomogene- 
ity of the polymers' distribution in space may play an 
important role. Postponing the consideration of such sit- 
uations to future work, we assume that the last term 
is negligible. In that case, Ily obeys a closed equa- 
tion and, noticing that II is proportional to the poly- 
mer contribution II to the stress tensor (where the co- 
efficient of proportionality includes the rigidity of a sin- 
gle molecule), we conclude that one can write a closed 
system of equations for riy and u (not including n). 
Furthermore, in the case n — const one can directly 
write Uij(x,t) = Bi(x,t)Bj(x,t), where B(x,t) is pro- 
portional to R(x,t). Using V • R — 0, one finds the 
following hydrodynamic equations for dilute polymer so- 
lutions 

d t u + (u ■ V)u = -Vp + vV 2 u + B ■ VB, V • u = 0, 
d t B + {u-V)B = (B-V)u-— , V-B = 0. (64) 

T 

The above equations were introduced in , based on the 
purely macroscopic consideration that the hydrodynamic 
equations admit the single-axis ansatz for Ily. The syn- 
chronization provides the microscopic basis, explaining 
the ansatz. It should be noted that, while the synchro- 
nization implies the single-axis ansatz, the converse is not 
true: one may have significant fluctuations in the end-to- 
end distances of individual polymers, non-contradicting 
the single-axis ansatz for the molecule-averaged quantity 

IV 

Thus, synchronization clarifies the derivation of hydro- 
dynamic equations for dilute polymer solutions, and can 
be used for a better understanding of the role of inhomo- 
geneity in the polymer distribution in space. 



XII. CONCLUSION 

In this work we derived new results on the behavior 
of small elastic objects in a turbulent flow. These results 
concern both the single particle behavior, and the behav- 
ior of many particles, even though it is always assumed 
that there is no interaction between the particles. The 
results can also be applied to chaotic flows satisfying the 
appropriate conditions of the decay of correlations. 

We first consider the behavior of small threads. The 
threads do not resist the average stretching by the flow, 
unless they reach the size comparable with their full 
length. Once this happens, the thread's size "stabilizes" 
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At < 1 At > 1 



FIG. 4: Sketch of the domain structure below (Ar < 1) and 
above (Ar > 1) the coil-stretch transition for polymers in a 
turbulent flow. The characteristic size of the ordered domains, 
which spontaneously arise below the coil-stretch transition, 
is the Kolmogorov length. Above the transition, if polymer 
back reaction is important, an even larger correlation length 
is possible [HEi]. 



at about the full length. The thread becomes similar to a 
rod, and it rotates in the same way as the vector connect- 
ing two fluid particles in the same flow. This continues 
unless the thread arrives at a region of the flow, where the 
local velocity gradient leads to a shrinking of the thread. 
The probability for the threads to shrink to a size con- 
siderably smaller than the full length is significant: it is 
known to decay only as a power law of the thread's ex- 
tension. We show that the exponent of this power law 
is universal and is given by d — 1, where d = 2, 3 is the 
space dimension. Notably, the answer is the same for any 
statistics of the velocity field, for which the Lagrangian 
velocity gradient has a finite correlation time. 

One could think that this universality has a simple ex- 
planation: the answer would have the same form, if the 
second end of the thread would perform a simple diffu- 
sion in a ball, centered at one end of the thread, with a 
radius given by the thread's length, and with a reflecting 
inner surface. However, the physical motion is quite dif- 
ferent: it is the logarithm of the distance that is driven 
by the noise and not the distance itself. The noise has 
a non-zero average, in contrast to diffusion. Finally, the 
precise ratio between the average and the dispersion of 
the noise, such as the one of the Kraichnan model, was 
shown to be important to get the correct answer. This 
becomes clear when one passes to consider polymers, that 
resist stretching by the flow with harmonic forces for sizes 
much smaller than their full length. Here too, a power 
law is known to describe the probability of the polymer 
shrinking to a size much smaller than the full length (it 
should be mentioned that this power law has a cutoff at 
the characteristic equilibrium size R co a - the coil's radius, 
where the thermal noise becomes important). However, 
this power law is known to be non-universal and to de- 
pend on the velocity statistics (and the polymer's relax- 
ation time). Nonetheless, we show that the exponent of 



the power law has a more universal form in the limit well 
above the coil-stretch transition, where only one constant 
Xd completely characterizes the velocity statistics. 

Our findings arc in accord with the recent experimen- 
tal findings in @ (see also [26j]), where a universal expo- 
nent was obtained for the power-law of the PDF of the 
polymer size in a shear flow. The shear flow is different 
from the turbulent flow we study here, where small-scale 
isotropy holds, and the generalization of our results to 
anisotropic situations such as shear, is postponed for fu- 
ture work. 

Furthermore, we show that different particles (poly- 
mers or threads) get synchronized by the flow that 
stretches them. Our conclusion concerning the synchro- 
nization of polymers has particular importance in view 
of the physical relevance for polymer solution. Below the 
coil-stretch transition, the polymers are coiled in most 
of the space, and the properties of the solution are very 
close to those of the solvent. The configurations of differ- 
ent coiled-up polymers get correlated because the same 
flow tries to stretch them, but the correlation is negligibly 
small. However, there are also regions where polymers 
are stretched to R > Rq, where Ro is some number much 
greater than R co a- The probability of such regions de- 
cays only as a power of R co u / Ro, and so these regions oc- 
cur in space rather frequently. The end-to-end distances 
in such regions are created by velocity-gradient fluctu- 
ations, for which the thermal noise is negligible during 
time-scales much larger than those needed for synchro- 
nization. Thus, in these regions the end-to-end distances 
of different polymer molecules are nearly identical (of 
course, the size of these regions is much smaller than the 
Kolmogorov length) . As the Lyapunov exponent of the 
flow grows, the " stretched regions" become more frequent 
in space, and above the transition they fill almost all the 
space. See the sketch of such a situation in Fig. |4j Thus, 
above the transition, there is a well-defined macroscopic 
field R(t, r) that varies over the scale of smoothness of 
the flow, and to which we also refer as domains. This 
field was originally introduced in Ref. [l4{ as a conse- 
quence of the hydrodynamic equations. Here, we give a 
microscopic meaning and interpretation to that field and 
discuss the implications for the hydrodynamic equations. 
The above predictions can be tested experimentally. 

We assumed in our analysis that the flow velocity is 
independent of the configurations of the polymers, so the 
stretching of the polymers is arrested by nonlinear elas- 
ticity. However, the most interesting case occurs when 
the stretching is arrested by the polymers' back reac- 
tion on the flow. Then there is a size Rback < Rmax 
such that, when polymers are extended to that size, 
their back reaction on the flow becomes significant. For- 
mally, the back reaction is realized via the polymer con- 
tribution to the momentum stress tensor. For example, 
if the polymer concentration is sufficiently large, then 
Rback <C Rmax, and one can apply the approximation of 
harmonic forces, where the contribution is proportional 
to RiRj, see [ll, . One can see that our derivation 
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of the power law still applies at R <C mm[Rb a ck, Rmax], 
where R decouples from a by the definition of Rbacki 
and our considerations can be repeated. As to synchro- 
nization, it is likely to hold without changes, but detailed 
study is needed. The first question would be, whether the 
orientations of different polymers still satisfy the same 
equation as the orientations of the end-to-end vectors of 
two fluid particles. If so, then - assuming there are no 
anomalies and the Jacobi matrix behaves in the usual 
way characteristic for chaos that was described above 
- the orientations will get equal, and one expects syn- 
chronization to hold. This conclusion is immediate if 
one could model the effect of back reaction on polymers 



by an effective radial resistive force that becomes impor- 
tant at R ~ Rback- The detailed study of whether our 
conclusions can be transferred to the important case of 
flows with back reaction, that includes the drag reduction 
problem, is the subject of future work. 
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